Experimental Demonstration of Squeezed State Quantum Averaging 
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We propose and experimentally demonstrate a universal quantum averaging process implementing 
the harmonic mean of quadrature variances. The harmonic mean protocol can be used to efficiently 
stabilize a set of fragile squeezed light sources with statistically fluctuating noise levels. The averaged 
variances are prepared probabilistically by means of linear optical interference and measurement in- 
duced conditioning. We verify that the implemented harmonic mean outperforms the standard 
arithmetic mean strategy. The effect of quantum averaging is experimentally tested both for un- 
corrected and partially correlated noise sources with sub-Poissonian shot noise or super-Poissonian 
shot noise characteristics. 
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The mean of a set of statistically varying, real and non- 
negative numbers, x = {xi,X2, ...x n }, is in general [1 : 




where r is an integer. The most commonly used means 
are the arithmetic mean and the harmonic mean corre- 
sponding to r = 1 and r = — 1, respectively [2]. These 
two kinds of means occur in several physical problems, 
and the actual mean being used to describe the physical 
system depends on the physical setup. E.g. the total re- 
sistance, R, of an electrical circuit consisting of n serially 
or parallelly connected resistors (with resistances Ri) is 
known to follow the arithmetic, R = Yl7=i anc ^ nar " 
monic, = Yli=i means, respectively (apart from 
a multiplicative constant). Similar mean laws can also 
be deduced for the total stiffness of a system comprising 
springs connected in series or in parallel. Likewise, one 
finds examples of the arithmetic mean and the harmonic 
mean in geometrical optics as well as in astronomy [3j H] . 

All these examples of the arithmetic and harmonic 
means are based on classical systems. In this Letter, we 
explore an example of the arithmetic and the harmonic 
mean in a quantum optical system. More specifically, we 
propose and experimentally demonstrate the arithmetic 
and harmonic means of the quadratures variances of dif- 
ferent quantum states using an optical system that is 
based on simple linear optics and homodyne detection. 
We average an ample supply of quantum states that ex- 
hibit either sub-Poissonian shot noise or super-Poissonian 
shot noise behavior, and we investigate the averaging pro- 
cedure for completely independent as well as partially 
noise correlated quantum states. Such averaging proto- 
cols allows us to stabilize the degree of squeezing (vari- 
ances) of n independent, fragile and possibly unstable 
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squeezed state resources. Therefore, besides being of fun- 
damental interest, such a protocol will find applications 
in quantum information and quantum metrology where 
stable resources of squeezed light are required [5 , 6 . We 
find that in terms of stabilizing the squeezed state vari- 
ance from a fluctuating set of different noisy squeezed 
light resources, the harmonic mean outperforms the stan- 
dard arithmetic mean. 




(n-l) Detections 



FIG. 1: Schematic illustrations of two arithmetic mean pro- 
tocols and the harmonic mean protocol (see text for details). 
SP: Signal processing; BS: Beam splitters; SQZ: Squeezing 
source; T(M): Trigger (modulator). Experimentally we ex- 
plore setup b) and c). 



Consider n independent quantum resources described 
by the quadrature variances Vi = (xf) — (xi) 2 where 
Xi, (i = l,2...n) are the amplitude quadratures. The 
goal is to construct the arithmetic and harmonic means 
of the variances without using any additional squeezed 
state resources. The arithmetic mean can be formed in 
two different ways as illustrated in Fig. [Ij a) By ran- 
domly selecting one of the n states from the resources or 
b) by interfering the n states on (n — l)-beam splitter [7 
and subsequently rejecting all outputs except one. Any 
one of these approaches yields the arithmetic mean for 
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the total variance: 

1 n 

Va = -5> (2) 

2=1 

Note that such a strategy is similar to the demonstration 
of universal quantum purification of continuous variable 
quantum information based on an ensemble of identical 
states |S]. 

The harmonic mean of the variances can be formed by 
using the setup shown in Fig.[T]3. Here the n states inter- 
fere on (n — l)-beam splitter, the amplitude quadratures 
of the (n — 1) outputs are measured and the results are 
used to drive the state. The state can be driven in two 
different ways: If the amplitude quadrature variances of 
the input states are known, the protocol is determinis- 
tic and the state is linearly displaced with an amount 
determined by the a priori variance information as well 
as the measurement outcomes. However, in a more rele- 
vant scenario where the variances are unknown, the state 
is probabilistically heralded based on the measurement 
outcomes: If the outcomes are arbitrarily close to zero, 
the state is kept, otherwise it is discarded. Importantly, 
the latter protocol is universal as it is independent on the 
input state. Either methods yield the harmonic mean for 
the variances: 

- = 1 E 1 - (3) 

1 = 1 

We note that a quantum optical version of the " resistor" - 
type harmonic mean for the amplitude quadratures, 
l/Vzr = J27=i V^n can be also implemented exactly 
solely using Gaussian operations, by replacing the array 
of beam splitters in Fig(T]: with quantum non-demolition 
interactions. This, however, requires the use of additional 
squeezing resources and is therefore not considered fur- 
ther in this paper. 
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FIG. 2: Comparison between the two different means. Vari- 
ances of the arithmetic and harmonic means a) for an in- 
creasing number of resources (one with variance V = 4 and 
V = 0.25 for the rest), and b) for a supply of five resources 
with an increasing number of noisy resources (with V = 4) 
and the rest being quiet with V = 0.25. 

By using the two averaging operations, it is possi- 
ble to stabilize the variances of n independent unstable 
squeezed state resources. As an example, let us consider 




FIG. 3: Schematic of the laboratory setup for the implemen- 
tation of quantum averaging. Squeezed states are prepared 
in optical parametric oscillators (OPOs) and additional noise 
can be added in amplitude modulators (AMs). The measure- 
ment induced harmonic mean operation and the verification 
are performed with homodyne detectors (HD1 and HD2). V\ 
and Vi are the input variances, Va is the arithmetic mean and 
Vh is the harmonic mean. BBS: Balanced beams splitter; SP: 
Signal processing; LO: Local oscillator. 



four noise resources with V = 0.25 and a single broken 
source with V = 4 (both variances are normalised to the 
variance of a vacuum state). The arithmetic mean trans- 
formation produces a single source with Va = 1 whereas 
the harmonic mean transformation produces a source 
with Vh = 0.31. Interestingly, the arithmetic mean pro- 
tocol stabilizes the sources to the shot noise level and the 
harmonic mean stabilizes it below the shot noise level. 
Moreover, if another source is broken (also with V = 4), 
then one gets an arithmetic mean of Va = 1.75 and a 
harmonic mean of Vh = 0.40, which means that the har- 
monic mean method is much less sensitive to the number 
of broken resources. The two means are compared in 
Fig. [2^i for different numbers of resources and a single 
noisy resource with V = 4, and in Fig. [2}d for different 
numbers of noisy resources with a total of 5 resources. 

We note that a stabilization procedure for an am- 
ple supply of continuous variable quantum informational 
states as well as qubit transformations were addressed 
theoretically in ref. [9] and refs. [TQl [11], respectively. 
Unlike these proposals, in our work we stabilize off-line 
resources without disturbing the actual quantum proces- 
sor. Furthermore we note that in contrast to Gaussian 
squeezed state distillation (which is not possible using lin- 
ear optics [H]), quantum averaging can be implemented 
with solely linear optical elements and feedforward. In 
addition, we stress that the squeezed state distillation of 
Non-Gaussian noise in refs. p~3j [14] was state dependent 
whereas our protocol is universal. 

We experimentally implement the two averaging trans- 
formations for a supply of two quantum states (n = 2) 
that exhibit either sub-Poissonian shot noise or super- 
Poissonian shot noise behavior. The schematics of our 
experiment is depicted in Fig. [3] Basically, the setup 
comprises three parts; a quantum state preparation part, 
a quantum protocol part and a verification part. We pre- 
pare two Gaussian states using optical parametric oscil- 
lators (OPO) followed by amplitude modulators (AM). 
The OPOs are bow-tie shaped cavities with type I peri- 
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odically poled KTP crystals [15j [16] . To produce a bright 
excitation of the squeezed beam and to enable a cavity 
phase lock, we inject two auxiliary beams into the cavity. 
Squeezed states are produced at 1064 nm by pumping the 
parametric process with a mode-matched beam at the 
wavelength of 532 nm. The relative phases between the 
pump and seed beams are locked to de-amplification in 
order to generate amplitude squeezed beams. To produce 
a whole range of different noise variances, the amplitude 
squeezed states are sent through amplitude modulators 
(AM) which are driven by electronic noise sources with 
variable modulation depths. Using such a combination 
of OPO and AM, the noise of the amplitude quadra- 
tures can be tuned from sub-Poissonian shot noise to 
super-Poissonian shot noise. To execute the protocol, the 
resulting beams interfere (with a visibility higher than 
98%) on a balanced beam-splitter (BBS). The relative 
phase is set such that the two beams add in quadra- 
tures. The two output beams of the beam-splitter are 
then measured by homodyne detectors (HD) with slowly 
varying local oscillator phases. The measurements are 
performed at the sideband frequency of 4 MHz with a 
bandwidth of 300 kHz, and the output signals were am- 
plified and digitized at 5 x 10 6 samples per second. Data 
bins associated with the measurements of the amplitude 
quadratures contains approximately 6 x 10 4 data points. 

The arithmetic mean is produced directly after the 
beam splitter by discarding the outcomes of HD1. To 
implement the harmonic mean protocol, we select the 
outcomes of HD2 based in the outcomes in HD1 with 
varying threshold values: If the data point of HD1 is 
lower than the threshold value, the corresponding data 
point of HD2 is kept, otherwise it is discarded. The 
state heralding process could in principle also be imple- 
mented electro-optically to generate a freely propagating 
averaged quantum state. However, to avoid such compli- 
cations, our conditioning is based on digital data post- 
selection. 

We present the experimental results for the arithmetic 
and harmonic mean protocols in Fig. [4j Fig. |4^i shows 
the results for two amplitude squeezed beams with vari- 
ances Vi = 0.64 ± 0.01 and V 2 = 0.90 ± 0.02. We 
see that for a success probability (ratio of the data 
kept after post selection to the initial data) of around 
0.10, the harmonic- mean method produces a state with 
Vh — 0.74 ± 0.02 whereas the arithmetic-mean method 
gives Va = 0.77 ± 0.02. In this particular case the im- 
provement of the amplitude noise using the harmonic 
mean method is very small and is basically within the 
measurement uncertainty. However, the superior perfor- 
mance of the harmonic mean with respect to the arith- 
metic mean is clearly manifested by considering the in- 
put variances of Vi = 0.62 ± 0.01 and V 2 = 1.83 ± 0.04. 
As shown in Fig. (4}d, in this case the arithmetic-mean 
method produces a state with Va = 1.22 ± 0.03 whereas 
the harmonic mean method produces a squeezed state 
for success probabilities lower than 0.60. With a success 
probability around 0.10, the harmonic- mean method gen- 
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FIG. 4: The arithmetic mean and harmonic mean as a func- 
tion of the success rate Ps for for two different supplies of 
resources. In a) the two input states are squeezed with vari- 
ances Vi = 0.64±0.01 and V 2 = 0.90±0.02, and in b) only one 
of the resources are squeezed with variance Vi = 0.62 ± 0.01 
while the other resource has a variance V2 = 1.83 ±0.02 above 
the quantum noise level (QNL). Theoretical predictions are 
represented by solid (harmonic mean) and dashed (arithmetic 
mean) curves. 



erates a state with Vh ~ 0.90 ± 0.03. This shows that 
the universal harmonic mean strategy can stabilize very 
fragile and unstable quantum noise sources below the 
quantum noise limit (QNL) against a source that sud- 
denly generates a large amount of classical noise. The 
solid curves in Fig. [4] represents the theory taking into 
account various imperfections of the setup. The error 
bars of the variances depend on the measurement error, 
which mainly are associated with the stability of the sys- 
tem over time. During a complete measurement (lasting 
3-4 hours) the standard deviation of the quantum noise 
level was 0.02 shot noise units. And the statistical errors 
which are due to the finite measurement time and the 
post selection process [T7] . 

In the above analysis and experiment, the n sources are 
completely independent and thus the amplitude quadra- 
tures are uncorrelated, C = (X\....X n ) = 0. We now 
consider the situation of partially correlated sources. Al- 
though the theory can be easily conducted for an arbi- 
trary number of sources, for simplicity we consider only 
the case of two sources (n = 2) with a quadrature cor- 
relation described by the coefficient C = (X1X2). Using 
the same setups as above, we find that the arithmetic 
mean is modified to Va c — Vl ~^ V2 C. whereas the proba- 
bilistic harmonic mean is approaching Vhc = 2 y^y~+ 2 c 
for a very narrow post selection interval. Through sim- 
ple inspection, we see that the effect of the correlations 
(C > 0) is a reduction of the mean values. Moreover, we 
easily find that Vh — Vhc if the correlations are unbi- 
ased and otherwise Vh > Vh c - We also stress that these 
protocols are, as above, universal in the sense that no 
a priori information about the input state is required to 
execute the transformations. On the other hand, if the 
noise is maximally correlated and unbiased (and this is 
a priori known), the noise can be removed perfectly and 
deterministically as shown in ref. [18 . 

We now test the universal averaging protocols exper- 
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FIG. 5: The influence of correlations between the resource 
states. In a) we plot the variances as a function of the success 
probability for different degrees of correlation and the results 
are summarized in b) for a success probability of 10%. The 
variances of the input states are V\ = 1.95 =b 0.04 and Vb = 
3.72 ± 0.07 for all realizations. Theoretical predictions are 
represented by solid (harmonic mean) and dashed (arithmetic 
mean) curves. 



imentally using input states with partially correlated 
noise. The two AMs were used to impose correlated noise 
onto the two states by employing a joint electronic noise 
generator. Due to the correlations, the noise interfere at 
the balanced beam splitter either constructively to pro- 
duce a highly noisy output state or destructively to pro- 
duce an output state with reduced noise. The output 
with reduced noise then directly serves as the output of 
the arithmetic mean protocol. To execute the harmonic 
mean, the output with increased noise is measured, post 
selected with different threshold values and finally used 
to herald the quadrature data measured with the veri- 
fying detector (HD2). The experimental results of the 
arithmetic and harmonic means for three different corre- 
lation coefficients are shown in Fig. [5^. For all implemen- 
tations the correlations were biased with the individual 
variances being Vi = 1.95 ± 0.04 and V 2 = 3.72 ± 0.07. 
We clearly see that as the correlation becomes stronger, 
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both the arithmetic and harmonic mean variances are re- 
duced. This trend is further illustrated in Fig. (5}d where 
the success probability for the harmonic mean is set to 
0.10. In all figures we insert the theoretical predictions 
(curves) based on the experimental parameters. This 
clearly demonstrates that the quantum averaging process 
is improved if the noises are correlated. 

In summary, we have extended the notion of arithmetic 
mean and, in particular, harmonic mean to the field of 
quantum optics. Several schemes of implementing the 
arithmetic and harmonic means of quadrature variances 
of a quantum state of light have been devised. Exper- 
imentally, we have demonstrated a probabilistic scheme 
for the implementation of the harmonic mean based on 
a measurement-induced operation, and the results have 
been compared to the results of a trivial arithmetic mean 
protocol. We found that the harmonic mean protocol is 
the best transformation for stabilizing squeezed state re- 
sources. 

It is interesting to note that the stabilization of 
squeezed states of light using the harmonic mean law can 
be readily extended to other media, such as squeezing of 
the collective spin of an atomic ensemble [19 , squeez- 
ing in Bose-Einstein condensates [20 and squeezing in 
plasmonic systems [15]. As the squeezed state has been 
shown to be the basic, irreducible off-line resource for 
universal state preparation [2TJ [22] and universal quan- 
tum computation with continuous variables [23j [24] , we 
foresee that the harmonic mean protocol might play a 
central role in future quantum informational and metro- 
logical technologies. 
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